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Abstract 

We perform a new, detailed calculation of the hadronic contributions to 
the running electromagnetic coupling, a, defined on the Z particle (91 GeV). We 
find for the hadronic contribution, including radiative corrections, 

10 5 x Z\ had .a(Mz) = 2740 ± 12, 

or, excluding the top quark contribution, 

10 5 x A had . a (5) (Mz) = 2747 ± 12. 

Adding the pure QED corrections we get a value for the running electro- 
magnetic coupling of 

1 



OQ.E.D.(M z 



128.965 ±0.017' 



(*) Work supported in part by CICYT, Spain 
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1. Introduction 

In a recent paper/ 1 ! hereafter to be referred to as TY-I, we have evaluated the hadronic contributions 
to the anomalous magnetic moment of the muon; specifically, a very precise determination of the piece 
involving the photon vacuum polarization function was given there. With a simple change of integration 
kernel (see below) this analysis can be extended to evaluate the hadronic contribution to the QED running 
coupling, oq.e.d.^)? m particular for t = M§; an important quantity that enters into precision evaluations 
of electroweak observables. We will find that we can produce a substantial improvement over previous 
determinations due to our use of complete and correct analyticity and unitarityl 2 ! properties (at low energy), 
and the high quality of recent Novosibirsk, LEP, and Beijing data. 
The running coupling constant may be written as 

oq.e.d.(*) = 6 ^ : , n(t) = e 2 iT rcn (t) (1.1) 

i + n(t) 

where e is the electron charge and 7T(i) rcn is the 1PI (one particle irreducible) vacuum polarization function, 
renormalized at t = 0. To lowest order we can write the shift in a as 

e 2 
aq.E.D.(t) = {1 + Aa(t)} — , Aa(t) = -e 2 77 ren (£). 

In fact, we will evaluate 77 ron (i) including the first radiative corrections, so the full (1.1) has to be used to find 
the effective coupling. However, we will follow current usage and will write, for the hadronic contributions 
77 h = e 2 iT. had 



rcn i 

A had a = -n h = -e 2 /7 h c a n d 



or, distinguishing between lowest order (index 0) and next order (index 1), 

Z\ (0) a = -7T (0) - _e 2 77had;(0) a(1) = _ rr(l) _ 2 ,jhad;(l) 

By using a dispersion relation one can write this hadronic contribution at energy squared t, 77/j(t), 
as 

-fl h (t) = -e 2 /7^(t) = -^ r As -^-, (1.2a) 

37T J 4m 2 S(S - t) 

with 

_, . a(e + e~ — > hadrons; s) (ri \, , , . Aira 2 , „, . 

R(s)= L, t _ qr^-4. ct() e+e-^/i + /i-; S =— — , 1.2b 

and the integral in (1.2a) has to be understood as a principal part integral. This is similar to the Brodsky- 
de Rafael expression for the hadronic contribution to the muon magnetic moment anomaly, 



o(h.v.p.) = / dsK(s)R(s), 



■Im'i 



K(s)=^k(s); K(s)= [ dx ■' J{ ~ 1 r) 



3ir 2 s ' J x 2 + (1 — x)s/m 2 L 

Therefore, we can carry over all the work from TY-I with the simple replacement 

K(s) — > ; r. 

For the coupling at the Z we will take t — M\. Because of this similarity with the g — 2 calculation, we will 
dispense with many discussions or details; they may be found in TY-I. Indeed, the present paper should be 
considered as a sequel to the former one. 

After the corresponding evaluations we find, to next to leading order in a, 

10 5 x Z\ had a(M|) = -10 5 x [77 n 0) (Mf ) + /^(Mf)] = 2740 ± 12, (1.3) 

or, excluding the top quark contribution, 

10 5 x Zi h adO (5) (M§) = 2747 ± 12. (1.4) 
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Adding the known pure QED corrections, the running QED coupling, in the momentum scheme is 

On E D (Ml) = . (1.5) 

q.e.d.v z) 128.965 ±0.017 v ' 

2. Contributions to the lowest order —II h in the energy range from threshold 

to 2 GeV 2 



2.1. The region s < 1.2 GeV 2 



r(0) 



To zero order in the cm. interactions we can write —11^ as a sum of contributions of different intermediate 
states in various energy slices. We start with the 2ir states for s < 1.2 GeV 2 . This we will subdivide in turn 
into two pieces: from threshold, 4m 2 , to 0.8 GeV 2 , and the higher energy piece. 

2.1.1. The region below 0.8 GeV 2 



We can express R^ ' in terms of the pion form factor, F^ 



R(°\s)= 1 -(l-^) 3/2 \FAs)\ 2 , (2-1) 



4 

where by m^ wc understand the charged pion mass. We can also relate F„ to the decay r + — ► D T ir + ir°. 
Consider the correlator 

n% v =i /Vze i ^(0|TV'+(a ; )K(0)|0) = (-p 2 ^ + P li p v ) H v (s) +p^n s (s), s=p 2 ; (2.2a) 

with V^ the weak vector current. Then neglecting isospin breaking (except for the phase space factor) we 
have at low s, 



v 1 (s) = 2irlmll v = ± 



1 (m^-m^f 



(m ir + + m^y 



3/2 



F*(*)\, (2.2b) 



and, on the other hand, v\ may be obtained from the experimental measurements of the decay r + — ► P T 7r + 7r°. 
To obtain F 7r (s) we will fit the recent Novosibirsk data' 3 ' on e + e~ — > 7r + 7r~ and the tau decay 
data of Aleph and Opal J 3 ' We will take into account, at least partially, isospin breaking effects by allowing 
different masses and widths for the p , p + resonances. Moreover, and to get a good grip in the low energy 
region where data are inexistent or very poor, we also fit F 7T (s) at spacelike sS 3 ' This is possible in our 
approach because we use an expression for F^ that takes fully into account its analyticity properties. ^ 2 ' To 
be precise, we use that the phase of F n (s) is equal to that of 7T7t scattering, in the elastic region, and then 
the Omnes-Muskhelishvili method. We write 

F n (s) = G(s)J(s). (2.3a) 

Here J is expressed in terms of the P-wave irir phase shift, S\, as 



«o 



J(s)=e i-^o)A !_M f 1 _±) cxp £/ d«'^i^- . (2.3b) 

So J \ s 0/ 7T 7 4 m2 S'{S' - S) 



Sl(s') 



s is the energy at which inelasticity starts becoming important (in practice, above the percent level); wc 
will take so = 1.1 GeV 2 in actual calculations. 
The exponential factor 



*»{-! as' s '^ 



S <s° 



^ Jiml s'(s'-s)j 

in Eq. (2.3b) guarantees that the phase of J(s) is equal to 6{(s) for s < so, hence equal also to the phase of 
F„. The rest is included so that J is smooth at s = So, and has the behaviour | J(s)\ ~ 1/s at large energies. 
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Because of this equality of the phase of J[s) and the phase of F w (s) below s = So, it follows that 
G(s) will be an analytic function also for 4to 2 < s < sq, so in the whole s plane except in a cut from s = sq 
to +00. If wc now make the conformal transformation 



r s^- V«o - s 



. , _ (2.4a) 

then, as a function of z, G will be analytic in the unit disk and we can thus write a convergent Taylor series 
for it. Incorporating the condition G(0) = 1, that follows from F n (0) = 1, and undoing the transformation, 
we have 



G(s) = l + Cl 



2\/«0 - V s -8 + 1 
_l,/sb~+ y/SQ- S 3 _ 


+ C 2 


( 5V^o- V s o - s 



(2.46) 



ci, C2, ... free parameters. Actually, only two terms will be necessary to fit the data. 

Next, to obtain J, and hence F w , we need a parameterization of (^(s). We can use the well-known 
effective range theory to write 



s 1 / 2 
cot*i(*)= 2fc3-(m 



s)tp(s), k = 



^s - Ami 



(2.5) 



where we have extracted the zero corresponding to the rho resonance. Now the effective range function ip(s) 
is analytic in the full s plane except for a cut for [—00, 0] and the inelastic cut [so, +00]. We can profit from 
this analyticity by making again a conformal transformation into the unit circle, which is now given by 



w 



Vso 



(2.6) 



V s + V s o - s 
We can therefore expand ip in a convergent series 1 of powers of w. Undoing the transformation wc then have 



5{(s) 



Arc cot 



2F (m " 



bo + h- 



Vso 



(2.7) 



We note that 60 = Const., 6^>i = would correspond to a pure Brcit Wigncr shape for the rho. By allowing 
for more terms in the expansion we are taking into account the known distortions of the Breit-Wigner shape 
due to the influence of the left and the inelastic cuts of tp. For the actual fits, only bo, b\, and m p are needed 
as parameters. 



The values of the parameters are obtained by fitting experimental data on e + e 



data on 



^ T 7T + 7r° decay, and data on F n (s) at spacelike s (ref. 3). We also include in the fit the value of the tttt 



P-wave scattering length, that we constrain at 



(38 ±3) x l(T 3 m^ 3 , 



(2.8) 



consistent with 7T7t scattering results as well as with current algebra calculations. For the free parameters of 
our fit wc find 



d = 0.23 ± 0.02, c 2 = -0.15 ± 0.03; b a = 1.062 ± 0.005, 61 = 0.25 ± 0.04; m p o = 772.6 ± 0.5 MeV . (2.9) 

We also find, as byproduct of our fit, the p° width as well as the mass and the width of the p + , the P-wave 
scattering length, and the mean square radius and second coefficient associated with the form factor of the 
pion: 

7> = 147.4 ±0.8, a\ = (41±2) x lO" 3 ™^ 3 , (r 2 ) = 0.435 ±0.002 fm, C?r = 3.60 ±0.03 GcV~ 4 (2.10a) 

and 

m p+ = 773.8 ± 0.6 MeV, T p+ = 147.3 ± 0.9 MeV . (2.10b) 

The x 2 /d.o.f. of the fit is 246/204 with only statistical errors, but improves to 214/204 when experimental 
systematic errors are included. 



It is to be noted that this series, as well as that in terms of 2 above, are quickly convergent in the region of interest 
for us here, which is mapped in segments contained in [—0.57, 0.24] inside the unit circles; see TY-I for details. 
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We have not included in this fit the experimental irir phase shifts (except for the scattering length) , 
as they are known to suffer from uncertainties associated with the method of extraction: 7T7t scattering 
cannot be measured directly. However, we have checked that adding them would not alter substantially our 
fit or parameters. Details of this, and other aspects of the calculation, may be found in TY-I, where also the 
results of separate fits to e + e" and tau decay data are presented. 

With the above parameterization of F^ we can evaluate immediately the corresponding contribution 
to 77. We find, with self-explanatory notation, 

-10 5 x 77^ 0) (M|; 2tt; s < 0.8 GeV 2 ) = 307.6 ± 2.2 ± 2.9. (2.11a) 

The first error is statistical, the second a combination of systematic (taking into account the correlations 
among the various sets of experimental data), and theoretical ones. 2 (2.11a) includes the (small) effect of 
lu — p mixing, evaluated with the standard Gounnaris-Sakurai method as in TY-I. 

Errors included in this work arc divided into statistical and systematic. Evaluation of the statistical 
errors is standard: the fit procedure (using the program MINUIT) provides the full error (correlation) 
matrix at the x 2 minimum. This matrix is used when calculating the corresponding integral for 77, therefore 
incorporating automatically all the correlations among the various fit parameters. 

In addition, for every energy region, we have considered the errors that stem from experimental 
systematics, as well as those originating from deficiencies of the theoretical analysis. The experimental 
systematics covers the errors given by the individual experiments included in the fits. Also, when conflicting 
sets of data exist, the calculation has been repeated, and the given systematic error bar enlarged to encompass 
all the possibilities. In general, errors (considered as uncorrelated) have been added in quadrature. The 
exceptions are explicitly discussed along the text. 

2.1.2. The 7T7T contribution in the region 0.8 < s < 1.2 GeV 2 

For the contribution in the region 0.8 < s < 1.2 GeV we integrate numerically the experimental data," 
and get 

-10 5 x 77^ 0) (M|; 2tt; 0.8 < s < 1.2 GeV 2 ) = 27.3 ± 0.3 ± 0.5. (2.11b) 

With the result above, 

-10 5 x 77^ 0) (M|; 2tt; 4m 2 < s < 1.2 GeV 2 ) = 334.9 ±4.1 (2.12) 

where both systematic errors (related to the same normalization uncertainty) have been added coherently. 

2.1.2. The 3tt, 2K, and other contributions in the region s < 1.2 GeV 2 

For the 3ir contribution we fit experimental data," with Breit-Wigner formulas (including the correct 
threshold behaviour) for the u), <f> resonances, plus a constant. We have two sets of experimental data; the 
difference between the evaluations with each of them is included into the systematic error. The x 2 /d.o.f. is 
63/60. The contribution to 77^ is, 

-10 5 x 77^ 0) (M|; 3tt; 9to 2 < s < 1.2 GeV 2 ) = 39.5 ± 0.3 ± 1.5. (2.13) 

The 27V states are treated in the same manner, fitting simultaneously e + e~ — ► K^Kg and e + e~ — > 
K + K~ data' 4 ' with the same Breit-Wigner parameters for the </>; the % 2 /d.o.f. is 84/82. For details we refer 
again to TY-I. We get 

-10 5 x 77<; 0) (M§; 2K; s < 1.2 GeV 2 ) = 41.6 ± 0.2 ± 1.3. (2.14) 

The contribution of 47r states is evaluated by numerical integration, with the trapezoid rule, of 
experimental data. <■" The systematic error includes the (estimated) difference between evaluations based on 
different sets of experimental data. This gives the result, 

-10 5 x 77f } (M|; 4tt; s < 1.2 GeV 2 ) = 2.6 ± 0.7 (2.15) 



If we had used data on e + e , but not on r decay, we would have obtained a slightly smaller number and a larger 
error: -10 5 x 77^ 0) (M|; 2tt; s < 0.8 GeV 2 ) = 306.5 ± 4.0 ± 4.3. We will take (2.11a) to be our best result here. 
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Finally, 5tt, 6ir, rjirw, . . . states contribute (0.3 ± 0.2) x 10~ 5 in this region. If we add all the 
contributions with s < 1.2 GeV we find, 

-10 5 x 7T^ 0) (M|; s < 1.2 GeV 2 ) = 418.9 ±4.6 (2.16) 

2.2. The energy range 1.2 GeV 2 < s < 2 GeV 2 

We have now a numerical evaluation obtained from a fit to inclusive e + e~ — ► hadrons experimental data:™ 
-10 5 x i7<; 0) (Mf ; 1.2 GeV 2 < s < 2 GeV 2 ) = 53.1 ± 5.3. (2.17) 

3. The lowest order -77^ in the energy range above 2 GeV 2 . The full 77^ (Mf) 

3.1. QCD calculations 

For the QCD calculations 3 we take the following approximation: away from quark thresholds, and for rif 
massless quark flavours with charges Qf, we write 



Jl(°)(t) = 3]Tq 2 j 1 + ^ + (1.986 - 0.115n/) {^j 



(3.1a) 



2 



- 6.64 - 1.20n/ - 0.005n) - 1.240 



(E f Qf)' 



To this one adds mass and nonperturbativc corrections. We take into account the quark mass effect for 
quarks with running mass fhi(s) which correct R^ ' by the amount, for each quark, 

-3Q 2 m 2 (s) | 6 + 28— + (294.8 - 12.3n/) ( — ) 1 s" 1 

+3Q 2 ^p- {-^ + f + (W - 10CO)) *■ j s~\ 

Finally, for the condensates we add 

0(l-i^)(a s G 2 )^Q 2 (3.1c) 



/ 

and 

247r 2 r a ~\ 

1-23-1 rm{il>iil>i). (3.1d) 

We neglect the condensates corresponding to heavy quarks (c, 6) and express those for u, d, s in terms of 
/ 2 m 2 , f^rrr^ using the well-known PCAC relations. The condensate contributions are negligible above 
s = 3GcV 2 . 

Eq. (3.1b) will be used when m 2 <C s. In practice this will mean that the contribution of the 
correction of order fhf/s 2 is less than 10 -5 . Near the threshold for heavy quarks c, b, t, i.e., when vf(s) <C 1 
(with Vi(s) = (1 — 4m 2 /s) 1 / 2 the velocity of the quark) we use a nonrelativistic QCD calculation (see refs. 7 
for details) in which the contribution of quark i is 

R?« = 3Q 2 [1 + 2 C0 ( S) ] 3 -ffl *°f ; (3.2a) 



See ref. 6 for the calculations of the various pieces. 
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a s {s) 



co(s) 



1 + 



(93 - 10n/)/36 + 7 e/3o/2 



a«(s), 



log- 



,1/2 



+o 4tt [ ° rriiCFOis 
To this we add the leading nonpcrturbative correction 

2ir(a s G 2 ) 



1-2 7E 



(3.2b) 



192mX 



and consider the effective threshold to occur when this overcomes the contribution (3.2a). 

In the intermediate region between vf <C 1 and m 2 <C s, we use the interpolation given by 
Schwingcrl 8 ! 



R 



Schw. 



ZQiVi(s 



3-vf(s) 



1 + C F 



TT 3 + Vi 


(* 


3\1 


1 


— 




3vi 4 


u 


4tt)\ 



(3.3) 



Note, however, that Schwinger's interpolation cannot be used for Vi — > as it underestimates Ri by a factor 
of 2. 

In the QCD calculations, the error labeled "Cond." is found by inserting the variation obtained 
setting quark and gluon condensates to zero, and that labeled A by varying the QCD parameter. Likewise, 
we label to, to the error obtained varying the mass m.j. If an error is not given it will mean that it falls 
below the 1CT 5 level. 

For the parameter A we take the recent determinations ™ that correspond to the value 

a s (M|) = 0.117 ±0.003; 
to be precise, we have taken (in MeV, and to four loops), 

A(s < to 2 ) = 373 ± 80; A(m 2 c < s < m 2 ) = 283 ± 50; A(m 2 b < s < m 2 ) = 199 ± 30; A(s > m t ) = 126 ± 20. 
For the gluon condensate we take (a s G 2 ) = 0.07 GeV . Finally, for the running quark masses we take 

m s (l GeV) = 0.188 GeV; m c {m c ) = 1.44 GeV; m b (m b ) = 4.3 GeV; m t {fh t ) = 174 GeV, 
and, for the pole masses, 

m c = 1.867 ±0.20 GeV; m b = 5.022 ± 0.060 GeV; m t = 174±5GeV. 
For the c, b masses, see refs. 10,11; for the t quark, ref. 12. 

3.2. The regions away from quark thresholds 

At the lowest energy region we find 

l5 



10° x 7T^ 0) (M|; 2 GeV 2 < s < 3 GeV 2 ) = 71.1 ± 0.5 (A) ± 0.4 (Cond.); 



(3.4) 



the justification of the applicability of QCD in this range is the agreement, within errors, of the QCD 
calculation with the e + e~ — > hadrons data, and with the more precise data coming from t — > v T + hadrons; 
this may be seen depicted in e.g. the plots of Aleph and Opal data in ref. 3 (more details may be found in 
TY-I). 

Apart from this region 2 GeV 2 < s < 3 GeV 2 , we can use the perturbative QCD formulas (3.1), 
(3.3) for the energy regions s > 3 GeV provided we stay away from heavy quark thresholds. We will thus 
get, excluding the J/ift, ip' resonances contributions (to be discussed below): 



-10 5 x 77^ 0) (M|; 3 GeV 2 < s < 3.7 2 GeV 2 ) = 259.1 ± 1.5 {A) 
(here the contribution of the error induced by the condensates is already negligible). Then, 
-10 5 x 77^ 0) (M|; 4.6 2 GeV 2 <s< 10.086 2 GeV 2 ) = 421.3 ±0.8 (A). 



(3.5a) 
(3.5b) 



(3.5c) 



(3.5d) 
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We will separate a region around M|, because we take the principal vale of the integral. We have 
thus, 

-10 5 x 77^ 0) (M|; 11.2 2 GeV 2 < s < 20 2 GeV 2 ) = 352.2 ±0.9; 

-10 5 x ZZf \M%; 20 2 GeV 2 < s < (M z - 3 GeV) 2 ) = 1668.9 ± 0.9; 
-10 5 x nl° } (M|; {M z - 3 GeV) 2 GeV 2 < s < (M z + 3 GeV) 2 ) = 29.2 ± 0.5; 

-10 5 x it^ 0) (M|; (M z + 3 GeV) 2 GeV 2 < s < 348 2 GeV 2 ) = - 794.5 ± 0.7. 
All the errors are due to the variation of the parameter A. Finally, 

-10 5 x 77^ 0) (M|; 360 2 GeV 2 < s < 400 2 GeV 2 ) = - 4.7 ± 0.3 

-10 5 x nl a) (M z ; 400 2 GeV 2 < s -► oo) = - 20.8 ±0.1. 

In particular, the total top quark contribution above threshold (360 2 GeV 2 < s — > oo) is -6.5. 

We note that part of the ranges contain some of the narrow resonances (tp, Y and T families). We 
will add their contributions individually later on. For the whole perturbative QCD contributions we have, 
adding Eq. (3.4) to Eq. (3.5), 

-10 5 x iT^ 0) (M|; 2 GeV 2 < s; pQCD) = 1982 ± 7. (3.6) 

Note that we have added the errors linearly as they stem from the same variation in the QCD parameter A. 

3.3. The thresholds regions 

We will make two types of calculations. In the first, we take experimental data (when possible, i.e., at the 
cc and bb thresholds); in the second, we take the contribution of the resonances lying below threshold from 
experiment, plus a background given by the contribution of the light quarks (evaluated with perturbative 
QCD, as above) and use nonrelativistic QCD to evaluate the contribution of the quarks whose threshold we 
are crossing. Of course, for the it threshold this is all we have. 

3.3.1. cc: J/ip, tp' and the continuum 3.7 2 GeV 2 < s < 4.6 2 GeV 2 

We split this into the contribution of the J/tp, ip' ', that we calculate in the n.w.a. (narrow width approxi- 
mation), and the rest. For the first we have, 

10~ 5 x (69.9 ± 4.5) [J ftp] 
10" 5 x (23.6 ±2.1) [V>'] 

For the remainder we have two possibilities: use a NRQCD calculation (see below) for the heavy 
quark, which gives 

10~ 5 x (73.2 ± 0.3 (A)) [uds]; (QCD; 3.7 2 GeV 2 < s < 4.6 2 GeV 2 ) 
10~ 5 x (66 ± 13 (m c )) [cc]. (NRQCD) 
Sum : 10~ 5 x (139 ± 13) 

Total: 10~ 5 x (233 ± 14) (QCD+NRQCD). 

Otherwise, we use experimental data above 3.7 2 GeV : 

10~ 5 x (111.8 ± 0.6 ± 5.5) (Exp., BES) : 3.7 2 GeV 2 < s < 4.6 2 GeV 2 
Total: 10~ 5 x (205.3 ± 7.4) (Exp., BES) 

(NRQCD) refers to the nonrelativistic QCD calculation with Eq. (3.2); sec TY-I and refs. 7 for 
details of this type of calculations. BES are the experimental data from ref. 13. The first error for them is 
the statistical, the second the systematic one. 

We give a few more details on the calculation with NRQCD, as it is the model for the other threshold 
regions. The method (QCD+NRQCD) consists in separating the u, d, s contribution; the cc one is then 

-7- 
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treated as follows. If a resonance is below the channel for open charm production, which is set at s = 4to 2 
(with c the pole mass of the c quark), then it is considered as a bound state, and treated in the n.w.a. 
Above cc threshold, one uses nonrelativistic QCD. For our choice of c quark mass, both J/ip and ip' should 
be considered to be below threshold. 

The reasonable agreement, within errors, between the (QCD+NRQCD) result for the cc contribu- 
tions in the region 3.7 2 GeV 2 < s < 4.6 2 GeV 2 , 1CP 5 x (139±13) and the result obtained using experimental 
data only, 1CP 5 x (112 ± 6), gives one confidence to use the same theoretical method of calculation for the 
other thresholds where the quality of the experimental data is poorer, or these data are lacking. However, 
for the cc region we consider that the results based on experimental data are the best and thus write 

-10 5 x 77^ 0) (M|; 3.7 2 GeV 2 < s < 4.6 2 GeV 2 ) = 205 ± 7. (3.7) 

3.3.2. bb: T, V and the continuum 10.086 2 GeV 2 < s < 11.2 2 GeV 2 

For the region around bb threshold we can repeat the calculations as above. First, we add the contribution 
of the resonances below bb threshold, that we calculate in the n.w.a.: 

1CT 5 x (5.8 ±0.2) [T] 

10~ 5 x (2.1 ±0.1) [T] 

For the continuum we find, for a b quark pole mass or 10 ' mf, = 5.022 ± 0.060 GeV, 

10~ 5 x (57.9 ± 0.1 {A)) [udsc]; (QCD; 10.086 2 GeV 2 < s < 11. 2 2 GeV 2 ) 
10" 5 x (8.7 ± 0.6 {m b )) [bb]. (NRQCD) 
Sum: 10~ 5 x (66.6 ± 0.6) 

If we had estimated the bb contribution saturating with the resonances Y", . . . T v , with electronic 
widths as given in rcf. 14 we would have got 5.2 ± 1.2 instead of the value 8.7 ± 0.6 that we found with the 
NRQCD calculation. We choose this last as our preferred value and write thus 

-10 5 x fl^ ] {M 2 z ; 10.086 2 GeV 2 < s < 11.2 2 GeV 2 ) = 75 ± 1. (3.8) 

3.3.3. It threshold: T bound states and the continuum 348 2 GeV 2 < s < 360 2 GeV 2 

The bound states produce a negligible contribution; for the ground state, a second order QCD calculation^ 11 ] 
gives r(T -+ e+e~) = 12. 5± 1.5kcV and thus the contribution to -10 5 x 77^ 0) is of -0.11. For the threshold 
region, a NRQCD calculation gives, for the t quark contribution, —0.47, while the udscb one is —1.41. (Note 
that in this calculation we are neglecting electroweak interactions, so we treat the t quark as if it was stable). 
All together, we find 

-10 5 x 7T^ 0) (M|;t thresh.) = -2. (3.9) 

The error is negligible. 

The total contribution of the threshold regions is thus 

-10 5 x 77^ 0) (M|; c, 6, t thrcsh's.) = 278 ± 7. (3.10) 

3.4. The lowest order i7^ 0) (Mf) 
Adding all the contributions to 77^ ' [M\) we get 

10 5 x Afl & u{Ml) = -10 5 x 77^ 0) (M|) = 2732 ± 12, (3.11) 

or, excluding the top quark contribution, 

10 5 x A^a {5) {Ml) = 2739 ± 12. (3.12) 
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4. The radiative corrections, —tl\j\ the full — tl^ +1 ' ; aQ.E.D.(-Wf) 

4.i. -np 

We have next the contribution of intermediate states containing a photon. At low energy (s < 1.2 GeV 2 ) 
we evaluate them individually, and at high energy (s > 1.2 GeV ) with the parton model. For the second 
we have a contribution equal to the zero order one (for which we take the result of the previous section) 
multiplied by the factor 

'YjfQf 3a 

This gives 

-10 5 x n£\M z ; s> 1.2 GeV 2 ) = 1.4 ±0.1, (4.1) 

the error depending on what one does in the quark thresholds, especially around the narrow resonances 
(J/7T, ip' , T, T'). For the low energy region we repeat, with obvious changes, the analysis of TY-I. Only the 
processes 7r + 7r~7, ir°j and 777 produce effects at the 1CP 5 level (respectively, 3.4 ±0.8, 2.9 ±0.2 and 0.8 ±0.1, 
in units of 10 -5 ). The first is evaluated in the narrow width approximation, or with a detailed calculation 
using theoretical formulas that relate 7T7T7 to 7T7t, and taking into account experimental cuts; the details may 
be found in TY-I. Both methods give essentially the same result. The other two are evaluated in the narrow 
width approximation, dominated by the p —^ 7r°7, lo — > 7r°7 and (j) — > 777 contributions. In addition, the low 
energy (s < 0.7 2 GeV 2 )7r°7 is calculated with a phcnomcnological coupling 7r°77, adjusted to reproduce the 
decay ir° — > 77. Again, the details are given in TY-I. 
Adding all of this to (4.1) we find 

-10 5 x 77^ 1) (M|) = 8.5±0.9. (4.2) 

We summarize our results in Table 1: 
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Channel 


Energy range 


Method of calculation 


Contribution to 

-10 s x Z\ha d a 


7T 7T _ 


s < 0.8 GeV 2 


Fit to e + e~ + r + spacel. data 


307.6 ±3.6 


7T + 7T~ 


0.8 < s < 1.2 GeV 2 


Fit to e + e _ data 


27.3 ±0.6 


3tt 


s < 1.2 GeV 2 


B.-W. + const, fit to e + e~ data 


39.5 ±1.5 


2K 


$ < 1.2 GeV 2 


B.-W. + const, fit to e + e~ data 


41.6 ±1.3 


4-7T, 57r, rjivn, . . . 


s < 1.2 GeV 2 


Fit to e + e~ data 


2.9 ±0.7 


Inclusive 


1.2 GeV 2 < s < 2 GeV 2 


Fit to e + e~ data 


53.1 ±5.3 


Inclusive; uds 


2 GeV 2 < s < 3.7 2 GeV 2 


Perturbative QCD 


330.2 ± 2.4 


J/4,, y/ 




n.w.a. 


93.5 ±5.0 


Inclusive 


3.7 2 GeV 2 < s < 4.6 2 GeV 2 


Fit to e + e~ data 


111.8 ±5.5 


Inclusive; udsc 


4.6 2 GeV 2 < s < 10.086 2 GeV 2 


Perturbative QCD 


421.3 ±0.8 


T, V 




n.w.a. 


7.9 ±0.2 


b quark thresh. 


10.086 2 GeV 2 < s < 11.2 2 GeV 2 


Pert. + Nonrelativistic QCD 


66.6 ±0.6 


Incl.; udscb(t) 


11. 2 2 GeV 2 < s < co (except t thr.) 


Perturbative QCD 


1230.3 ±3.4 


t quark thresh. 


348 2 GeV 2 < s < 360 2 GeV 2 


Pert. + Nonrelativistic QCD 


-2.0 ±0.1 


7+ hadrons 


Full range 


Various methods 


8.5 ±0.9 



Tabic 1 



Summary of contributions to Z\hadtv = — Ilh- "B.-W. + const." means a Breit-Wigncr fit, including the 
correct phase space factors, plus a constant; note that only for the four narrow resonances J/tp, ip 1 ] T, T' 
we use the n.w.a. The errors are uncorrelated except those for QCD calculations (that have to be added 
linearly) and those for the 2tt states (see text). For the details of the final states 7+ hadrons we refer to 
TY-I. 



4.2. The full ZihadCt; the QED coupling on the Z; discussion 

Adding then (4.2) to the lowest order expression we get the final result 



•10 5 x 



r(0) 



t(i)/ 



n^'(M'i) + n?> {Ml 



2740 ±12, 



10 5 x Z\ had a(M|) = 
or, excluding the top quark contribution, 

10 5 x Z\ had a (5) (M|) = 2747 ± 12. 
The pure QED corrections amount to (see, e.g., ref. 15) 

-10 5 x 77q ED (M|) = 3149.7687. 
Adding this to (4.3) and using (1.1) we get the value for the running electromagnetic coupling 



a Q . E .D.(M|) = 



1 



128.965 ±0.017 



(4.3) 
(4.4) 
(4.6) 

(4.7) 
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Whcn comparing with other determinations we will restrict ourselves to those performed after the 
data from Novosibirsk! 3,4 ] and Beijing! 13 ' have become available; these experiments increase the set of data by 
almost one order of magnitude, and have much better precision than the older ones. Thus, one may consider 
older determinations as superseded. So we compare our results with the determinations of refs. 15,16. We 
have, 

f 2743 ± 19 / 2765 ± 21, (MOR) 
10 5 x Zi had a (5) (^I) = I 2761 ± 36, (BP) (4.5) 

[2790 ±40, (J) 

In the MOR determination, the two values depend on the method of calculation used. As a general rule, 
comparing any two of the results quoted above, the largest reason for a reduced error bar is related to a 
wider use of perturbative QCD. For instance, the analysis of MOR uses perturbative QCD in the regions 
2.8 GeV 2 < s < 3.7 2 GeV 2 , and s > 5 2 GeV 2 . We have used pertubative QCD in the region 2 GeV 2 < s < 
3.7 2 GeV , justified in view of its agreement with the precise new experimental data (as discussed in the 
text), and the regions above the cc and bb thresholds. Here, the use of the calculations incorporating the 
exact effect of the quark masses has allowed to get a precise determination for 4.6 2 GeV 2 < s < 5 2 GeV 2 , 
where the experimental data are not very precise, but (again) perfectly consistent with QCD. 

In conclusion, we have performed a detailed evaluation of the hadronic contributions to the running 
electromagnetic coupling, obtaining a substantially reduced error bar. The ingredients are the following: 
First, we use Novosibirsk (e + e~) and LEP (r) data to fit the 2tt contribution. Invoking the analyticity and 
unitarity properties of the pion form factor allows to include spacclikc data also, improving the compatibility 
of the e + e~ data with the results from r decay, and reducing the corresponding error. Second, the low 
energy 37r and IK states have been considered individually, after the latest Novosibirsk data on the ui and 
(j> resonances. We perform a full-fledged fit, including the exact threshold factors. Third, we have used 
perturbative QCD in the region s > 2 GeV (away from quark thresholds). In particular, the recent LEP 
t — > v T + hadrons, and BES e + e~ — ► hadrons data justify the QCD result for 2 GeV 2 < s < 3 2 GeV 2 , 
implying the largest part of the error reduction. We also use the Beijing! 13 ! data, thus gaining precision, for 
the contribution in the energy range 3.7 2 GeV 2 to 4.6 2 GeV 2 . Last but not least, the next order radiative 
corrections have been taken into account. This is essential for our calculation as the radiative contribution 
is indeed of the same order of the final error bar. 
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